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Abstract 

We report on results of a theoretical study of the adsorption process of a single carbon oxide 
molecule on a Platinum (111) surface. A four-component relativistic density functional method was 
applied to account for a proper description of the strong relativistic effects. A limited number of 
atoms in the framework of a cluster approach is used to describe the surface. Different adsorption 
sites are investigated. We found that CO is preferably adsorbed at the top position. 
PACS: 31.10+z, 31.15.Ew, 31.15.Ne, 31.30.Jv 
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I. INTRODUCTION 



The interaction between atoms and surfaces is essential to understand adsorption pro- 
cesses. In order to determine accurate structures, vibration frequencies, electronic properties 
and forces ab initio calculations are indispensable. Density functional theory (DFT) nowa- 
days provides a powerful tool in nearly all branches of fundamental research in microscopic 
regimes (see e.g. ]]], 0, [3], U). It also facilitates a detailed theoretical investigation of prob- 
lems in surface physics and catalysis on a quantum-chemical level of theory ||, particularly 
the adsorption process of ad atoms on surfaces within the framework of a molecular orbital 
picture. Due to the adsorption process on the surface the translation symmetry parallel to 
the surface is violated. From a strict theoretical point of view methods of solid state physics 
are no longer applicable for this problem. In order to treat such systems approximately one 
can use two different approaches; one is the solid state method using super cells or a slab 
model, and the other is the molecular method using clusters. As the size of the super cell or 
cluster increases, both methods should converge and yield the correct experimental value. In 
this paper we follow the molecular approach in which the physical system is approximated 
by a big molecule or cluster. Results of ab initio calculations are comparable with experi- 
mental results for small molecular systems and/or clusters. Such a cluster model is of course 
only an approximation for the description of unlimited surfaces and solids. Nevertheless, as 
long as the physical process which we are going to describe is strongly localized, a cluster 
scheme may be an appropriate description. The adsorption process of a single atom or a 
simple molecule is considered to be such a local effect. 

One of the most extensively studied systems in this respect is the adsorption of carbon 
oxide (CO) on a Platinum surface. However, no good quantitative description has been 
achieved so far. The reason might be the lack of a proper description of the strong relativistic 
effects resulting from the 5c? band of Pt, which lead to a different behaviour in the presence 
of Pt than on the surfaces of the homologue elements Ni and Pd || . 

We have already tested this cluster approach on simpler systems successfully. The ad- 
sorption of a single Na and Ba atom at the Na(llO) and the Ba(llO) surface, respectively, 
have been studied 0. The results were reasonable and encouraged us to apply our method 
to the more complex system widely used in catalytic processes: the adsorption of CO on a 
Pt(lll) surface. 
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In this paper we present the results of four- component relativistic density functional 
calculations for the adsorption of CO on a Pt(lll) surface. Recently, our computer code 
has been developed to a new quality which allows to calculate fragmentation energies and 
geometries of bigger clusters with local and semi-local density functionals. Compared with 
scalar relativistic approaches, the spin-orbit coupling is fully included in our formalism. 



II. WHAT IS KNOWN FOR THE SYSTEM CO AT PT(lll) SO FAR? 

The adsorption of CO at a platinum surface seems to be one of the most analyzed ex- 
perimental systems. Apai et al. || showed that the CO molecule is adsorbed linearly with 
the C bonding towards the surface which can be explained by the fact that the HOMO 
(5a) and LUMO [2tt*) of the CO are located at the C atom. So this site is preferred for 
a binding with a surface. The molecule does not dissociate after adsorption and the main 
part of the binding is dominated by the d band of the platinum surface. The CO remains 
to be effectively neutral. 

A theoretical explanation of the binding formalism was given by Blyholder ||. He studied 
the binding of CO at metallic surfaces by using a Hiickel method. Within this framework 
the results show a formal charge transfer from the bonding 5a orbital of the CO towards the 
metal. On the other hand there is a charge back transfer from the metal towards the anti 
bonding 2ir* of the molecule. This exchange of charge explains the fact that the molecule is 
still neutral after adsorption. 

With this theory one is able to distinguish the adsorption site at the surface by the 
vibrational frequencies of the adsorbed molecule: Effectively there is a charge transfer in 
the molecule from the bonding 5a to the anti bonding 2tt* orbitals. The occupation of 
the anti bonding orbital leads to a reduced electronic charge between the C and O atom. 
This results in a weakening of the bonding, and therefore the vibrational frequency of the 
adsorbed molecule is lowered. The strength of this effect depends on the adsorption site. A 
molecule adsorbed at a hollow position has a stronger overlap between the metal band and 
the 27r* orbital than at a top position. This is due to the fact that the molecule is closer 
to the surface. Therefore the frequencies are lower at hollow positions and one is able to 
distinguish the adsorption sites || |10] . 



Ertl et al. [11] reported their results of the adsorption of CO at the Pt(lll) surface. 
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They stated that the strongest bonding arises at a hollow position which they erroneously 
explained by the knowledge of similar systems like Ni and Pd. With higher coverage the 
bridge position is occupied building a a/(3) x y(3)/i£30° overlay structure. In the regime 
of zero coverage they estimated the binding energy to be 1.43 ± 0.04 eV. In the following 
years a lot of experimental papers appeared, discussing the binding energy in this limit of 
low coverage. The preferred adsorption site is now assumed to be the top position followed 



by the bridge position |12], [L3|, [14], |I5|, [16], [T7|, [18], [19|| . The results show an adsorption energy 
of 1.3 to 1.6 eV. All these results have one thing in common: The adsorption energy is 
calculated via the modified Clausius-Clapeyron formula of thermodynamics: Pressure and 
temperature are varied whereas the coverage is kept constant which is monitored by the 
change of the work function. The resulting curve is a measure for the adsorption energy. 
The adsorbed molecules form a dipole layer on the surface resulting in a change of the work 
function which is proportional to the coverage. In the regime of low coverage this assumption 
may not be valid. 

Newer results of Yeo et al. [p_0| show a significantly stronger binding energy of 1.89 ± 
0.20 eV at low coverage. They used the new Single Crystal Adsorption Calorimetry (SCAC) 
method which allows the direct measurement of the IR-photon emitted during adsorption. 

Somorjai et al. were able to estimate the binding distances of the adsorbed molecules. 
They found a distance of 1.85 ±0.1 A between the C atom and the surface at the top position 
and a value of 1.55 ± 0.07 A at the bridge position. The bridge position was observed at a 
coverage of |. The binding distance within the CO molecule after adsorption is measured 
to be 1.15 ± 0.05 A [pi , p2]| . This is a small elongation of the intra-molecular distance. The 



free CO molecule has a value of 1.12 A |23 



On the theoretical side there are a lot of different approaches in the estimation of the 



binding energy. Most of them use a cluster approach U% BR E6L p7L pq, B9, ROl 31, 32 



33| with different methods like modified Hiickel methods or pseudopotentials and different 
numbers of atoms. Some [27, E_8], |2j], |30| used only 1-2 atoms for the simulation of the 
surface. The estimated binding energies at the top position are between 0.19 eV and 4.21 eV 
reflecting the fact that non-relativistic pseudo-potential or Hiickel models are not able to 
give reliable results for the adsorption energy. The main reason for this seems to be the lack 
of reproducing a good behavior of the d orbitals and therefore the d band of Pt. Models 
with less than 10 atoms in the cluster will probably only lead by chance to good results. The 
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Atom Superposition and Electron Derealization model from Ray and Anderson [|25[] based 
on extended Hiickel method combined with atomic one particle energies and atom-atom 
repulsion, on the other hand, leads to acceptable results with a binding energy of 1.66 eV 
for the top, 1.26 eV for the bridge and 1.11 eV for the fee hollow position followed by 1.03 
eV for the hep hollow position. Thus the fee hollow position is slightly preferred among the 
both possible hollow positions. 



Other authors used slab models |34|, |35], || or a super-cell model |37|. However, all these 



methods were not applied in the regime of low coverage. They work with an effective coverage 
between | and |. At these coverages the lateral interactions of the adsorbed molecules are 
not negligible. Yeo et al. [^(J showed a decrease of the energy by 0.2 — 0.4 eV within this 
regime. Thus the results they get are in good agreement with experimental results but are 
not valid at low coverage. 

Brako et al. studied the on-top chemisorption of CO molecules as a function of the 
coverage and the lateral relaxation of the platinum surface using a classical Hamiltonian 
with an analytical model potential |38 . 



III. THE THEORETICAL METHOD 

We all know that an exact solution of the relativistic many-particle Dirac equation is 
not possible. Large-scale multi-configuration or coupled cluster approaches are the best 
approximations for atoms so far. This high quality relativistic calculations are not yet 
possible for more than diatomic molecules. 

For the size of systems we are discussing here only relativistic pseudo-potential or density 
functional calculations are available. The code which we are using |39| has its first roots 



in the relativistic version of the DVM approach |]40|| . Although the general structure of the 
density functional code has been described in several publications we nevertheless repeat it 
here in short because a number of details have been modified. 

Following the Hohenberg-Kohn [|T1J and Kohn-Sham theorems [fJ2[] and starting from 



the no-pair approximation jfij and also neglecting the minor important contributions from 
spatial components of the four-current j v (r) = the total energy (in atomic units) 
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may be written as 



with the density 



E[p] =T S + E N [p] + E H [p] + E xc [p] (3.1) 



j°(f)=p(f)= ^Wfc(r) (3.2) 

— mC 2 <£k<SF 

obtained from N one-particle Kohn-Sham Dirac spinors. This leads to the corresponding 
relativistic Kohn-Sham equations 

{T s + V N (f) + V c (f) + V xc (f)} ^(f) = eiVi(r). (3.3) 

Here T s = cctiPi + (/3 — l)c 2 is the relativistic operator of the electronic kinetic energy, 
V (r) represents the nuclear potential assumed as a sum of point charges 

W) = £#^' (3-4) 

V xc (f) is the exchange-correlation potential derived from the exchange-correlation energy 
E^ 



J xc 



yxc^^SJ^M (35) 
6p[f) 

For the self-consistent solution of equation (|3.3|) we use the relativistic local density 



approximation [13, EH] (RLDA) together with the Vosko, Wilk and Nusair parameterization 



4_6| for correlation. Nonlocal corrections, with the relativistic form (RGGA) |^7[ of Becke's 



GGA approximation p8| for exchange and Perdew functional of correlation |49j (B88/P86)is 
used perturbatively. Optionally we also have a relativistic extension of Perdew and Wang 
5tJ (PW91) for the xc-energy. The Coulomb potential V c (r) is given by 



V c (f)= I -^Xd 3 r'. (3.6) 
J |y?_ r /| 

Numerical evaluation of the Coulomb matrix elements requires the evaluation of the Coulomb 
potential at the grid points, which leads to a large number (n 2 p; n: number of basis functions, 
p: number of grid points) of nuclear attraction type integrals. It is possible to facilitate the 
n 2 p problem by approximating the exact density with an expansion in one center fit functions 

m = EEEE d^ax^rr^Aa), (3.7) 



where f^iXa) is the radial density of the wave function r of atom a. Xa is a radial vector 
centered at each atom x~a = r — R a . d^ l m are the expansion coefficients of the approximate 
density which are obtained via a least square fit to the density 



(p(r)-p(r)) 2 d 3 r , (3.8) 
or via minimizing the Coulomb energy of the difference density []5Tj 

■d 3 r'd 3 r . (3.9) 



(p(f) - p{r)){p{r) - p(r)) , :! 



For simpler molecules, that means without complex density structure, it was sufficient to 
use only the monopole part of the atomic wave functions to build the approximate density. 
With the open d shell of Pt it is necessary to use higher moments in the expansion 
Actually we use values up to Z = 2. Using this procedure the least square value of eqn. 
was less than 5%. 

With the use of the approximate density ( |3.7|) we can reduce the three dimensional integral 
in equation (|3.6|) to one dimensional integrals: 

A M a L r I 
a r=l 1=0 m=—l 



with 



4:71 1 

T/rm/ n yr 



.0 Xa 

In order to solve the Kohn-Sham equations ( |3.3| ) we use the Molecular-Orbital Linear 
Combination of Atomic Orbitals (MO-LCAO)-method and expand the molecular orbital 
wave functions in symmetry-adapted wave function Xji which themselves are expanded in 
atomic orbitals Xn^i^) which are atomic four-component Dir ac-Spinor s |52|, [53] 



Mr) = Y^Xjtfcji (3.10) 

3 

Xjif) = ^Xn v (f)d nvj (3.11) 

n v 

with n u = (u,n, K,m). Here Xn^i^) are the relativistic four-component atomic basis func- 
tions centered at the atoms. Using this symmetry adapted atomic basis set implies the 



non pair approximation The symmetry coefficients d Uu j are calculated by the program 



TSYM of Meyer [p2| , |53| . The expansion coefficients Cij of the molecular wave functions 
together with the one-particle energies 6j are the results of the SCF calculation of equation 
(P-3|). Inserting ( |3.10| ) into ( |3.3| ) gives the matrix secular equation in the symmetry orbital 
basis x, 

HC = eSC (3.12) 

where H and S are the Hamiltonian and overlap matrices respectively, C is the coefficient- 
matrix (cij) and e is the eigenvalue- vector (e^). The matrix elements, Hij and Sij are 



evaluated numerically using the multi-center integration scheme by Baerends et al. |5J, |55 
In this scheme the whole molecular volume is spatially partitioned into so called Voronoi cells 
and each cell is divided into an atomic sphere around the atomic nucleus and polyhedra. This 
method was extended by a special transformation to account for the relativistic behavior of 
the wave function near the nuclei ||3^ . 



IV. THE FROZEN CORE APPROXIMATION 



In order to reduce the computational requirement one very popular way is to use either 
pseudopotentials or the Frozen Core approximation (FCA). Both approximations are based 
on the idea that chemical binding is realized via the outer valence orbitals only. In our 
calculations we used the Frozen Core Approximation. In FCA it is assumed that the inner 
orbitals do not change their form during the adsorption process so that they can be kept 
fixed in the self consistence procedure. This reduces the numerical effort for the heavier 
elements like Pt drastically. 

Since the FCA is a very popular approximation and used in almost all kinds of ab ini- 
tio calculations (from density functional to quantum chemical methods) there are many 
references in the literature |56], |57|, [58] for this method. 

Although the principal procedure of the FCA in all ab initio methods is the same, there 
are differences in practical implementations. Therefore we give a short summary of the 
FCA in our code. Compared with all electron calculations there are two differences in 
FC calculations: the orthogonalization of the valence orbitals to the core orbitals and the 
calculation of the total energy. 

The FCA is based on a decomposition of the Hilbert space into a core (c) and a valence 
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(v) part. Accordingly, the molecular wave functions can be written as: 

In order to ensure the constrained Ritz variation of the valence orbitals, a core-valence or- 
thogonalization has to be carried out. In our implementation we followed the idea of the 
normalization of the overlap matrix in the core-valence orthogonalization. The decomposi- 
tion of the Hilbert space results in the decomposition of the overlap and Fock matrices into 
core-core, core-valence and valence-valence parts which have to be orthogonalized to each 
other. 

In order to ensure the core- valence orthogonalization all contributions of core orbitals are 
removed from the valence orbitals. This is achieved by the explicit orthogonalization of the 
valence orbitals onto the core orbitals. 

\Xv) = \Xv) S VC S CC \xc) ■ 

Here S vc is the valence-core overlap matrix. As a result the transformed overlap and Fock 
matrices of the whole system are now of the form: 

v V ' 

— ^neglected 

We assume that the core orbitals of different atoms do not overlap. By the use of atomic 
orbitals the core orbitals of the same atom are by default orthogonal. Therefore the corre- 
sponding symmetry orbitals are also orthonormal and S cc = 1. 

Due to the chemical binding the relative change in the core potential is small. The change 
in the core wave functions which is a second order effect can therefore be neglected. This 
means that H' vc in eqn. ( |4.l| ) can be neglected. As a result of this approximation the secular 
equation (|3.12|) decouples into two equations, which are solved like all electron calculations. 
Both approximations define our implementation of FCA. 
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The total energy within this Frozen-Core approximation is now: 

A Ng A 



a=l ft=l a=l 

a=l J a=l J 

a=l fi^a a=l f3^a 

M i r ~ r 

+ E n"^-- p»V v c d 3 r- V xc p v d*r 



+E™+J2 ■ (4-2) 

a</3 

The symbols with a small v represent atomic potentials. The subindices c and v mean 
that the core and valence densities are used to calculate the potentials. The first two lines 
of equation ( |4.2|) consist of pure atomic contributions. As far as we do not change the 
basis functions these values can be neglected as long as we are only interested in binding 
energies and distances. In this case these core contributions cancel completely out. The 
third line describes only the core-core interactions of different atoms. These values can be 
approximated by the use of point charges: 



E / VaP?d 3 r + ±J2 [ ^p^ r = EE / 



A iv N a + vg) d*r 



A A rye I rv \ ry c \ 

a=l fl^a \ rL ct JX P\ 

The last two lines are exactly the same as in the total energy expression in an all electron 
calculation but only with the use of the new valence orbitals for the density. For V xc we 
still need the full density p = p c + p v . 



V. RESULTS 

Our aim in this paper is the ab initio calculation of the adsorption of CO at Pt(lll) 
surface. This process is one of the most frequently used catalytic processes but is not 
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understood in detail up to now. We have referred to the experimental findings and theoretical 
methods which were applied up to now in section II. 

As Pt is a very heavy metal where relativistic effects play an important role we at least 
have the appropriate method. Unfortunately the computational effort restricts the actual 
calculations to a very limited number of Pt atoms which can be used in a real calculation. 
Nevertheless we have achieved to converge clusters between 7 and 13 Pt atoms depending 
on the adsorption site of the CO molecule. To keep the calculation time in an acceptable 
range, we in addition used only a quasi minimal basis set which was optimized via the degree 
of ionization of the atomic wave functions used as basis functions. For Pt we use the Is 
to 6p wave functions with a degree of ionization of +0.2. As an effect this leads to valence 
functions which are slightly more localized than the neutral ones but are more suitable to 
reproduce the binding character of the bulk. For the wave functions with the main quantum 
numbers n=l to 4 in addition the FCA was used. 

A. The CO molecule 

For the calculation of the CO molecule a larger basis set was used. This system requires 
a very good basis to account for the strong chemical bonding and charge distribution within 
this molecule. 

The best results for the CO molecule were obtained with a minimal basis set consisting 
of the Is — 2p wave functions of the neutral atoms in combination with additional 3p, 3d 
wave functions with an degree of ionization of +4.0. In this molecule all electrons were used 
in the SCF process. 

Table | shows the results for the CO molecule, which are close to the experimental data. 
The negative value of the dipole moment indicates a major negative charge at the C atom. 
The binding energy is slightly too high, which may be attributed to an unpolarized treatment 
because current-dependent functionals are not available. 

B. The adsorption calculations 

In order to calculate the adsorption energy of CO on the Pt(lll) surface one first has 
to define the clusters which simulate the surface for the four possible adsorption sites of 



11 



the CO molecule. The internuclear axis of the CO molecule is assumed to be perpendicular 
to the surface with C towards the surface. The four possible adsorption sites are the top, 
bridge, fcc-hollow and hep-hollow of the (111) surface of Pt. The different adsorption sites 
are shown in Figures [TJ to As insert in the figures [I] to [| we present the Pt-clusters as 
molecular models in the view from above and in a perspective view for the four adsorption 
sites. The CO molecule can also be seen. The distances between the Pt atoms are kept fixed 
to their bulk values (do = 3.92 A) and the CO distance to the value calculated above. 

In the cases of top and the hollow positions the C^ v symmetry was utilized. For the 
bridge position only a C\h symmetry can be used. Figures [l] to [| show the results of the 
calculations. Because of the low symmetry of the bridge position a maximum of 7 atoms had 
to be used to get results in a realistic calculation time. The notation Pt m+n (m,n) denotes 
the number of atoms per layer: m atoms in the first layer and n in the second resulting 
inm + n total Pt atoms. The bulk was represented only by two layers. As shown in the 



literature, this number of layers seems to be necessary and nearly sufficient [6C, Km 32 



Results of our calculations for binding of the CO molecule on the four surfaces are given 
in Tables |Tl -[IV] for the GGA calculations according to Becke 1988 [4S] in combination with 
Perdew 1986 |^9[ (which in short term notation is B88/P86). 

The results show the experimentally observed binding order with the top position as the 
strongest one and therefore as the preferred adsorption site. Like Ray and Anderson |^S] we 
also find a slight preference of the fee position in the hollow positions. 

Within our method the binding energies are by about 10 — 20% too high compared with 
experimental results (see for instance the results of the CO molecule). Having this in mind 
our results of 2.3 eV binding energy at the top position compares well with the measured 
value of 1.86 ± 0.20 eV of Yeo et al. [^0f (see table U). A direct comparison of binding 
energies at other positions with experimental results is not available. The calculations here 
were performed with a single CO molecule only. In experiments an occupation of the other 
adsorption sites is only observed at much higher coverages which are not considered here. 

In Table |II| the calculated bond distances of the C atom (of the CO molecule) with 
respect to the surface at the energy minimum at the four possible adsorption sites are listed. 

In addition very informative is Table [TV]. Here the distances are calculated towards 
the next Pt atom of the surface for all four CO positions. For better comparison with the 
experimental results these values are converted to A. These results correlate with the binding 
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behaviour. The distance at top position is the shortest whereas the hollow positions have 
the greatest distances. A comparison with the experimental values shows a good agreement 
for the top position as well as the bridge position where our value is slightly shorter than the 
experimental value. One reason for this may lie in the effective coverage when the bridge 
position is experimentally observed. As already mentioned, the interactions of the adsorbed 
CO molecules are not negligible P^J within this coverage. 

Table [V| presents a comparison of our values for the top-position with other theoretical 
results. These values demonstrate that the bond energies still do scatter up to 1 eV. The 



same applies for the various theoretical results for the hollow-position in Table [VI . 

Table |VII| shows the electronic charge distribution of the CO molecule as a result of a 
Mulliken population analysis of our calculations. The charge of the oxygen atom remains 
more or less constant and independent of the adsorption process and adsorption site whereas 
the carbon atom changes its charge by small amounts. A more detailed analysis of the 
electronic charge distribution shows that the charge between the two atoms C and O is 
reduced in the direction off the molecular axis. 

It is further interesting to see that the CO molecule as a whole has almost no charge 
transfer at the top position which has the strongest binding whereas all three other positions 
show a negative charge transfer from the cluster to the molecule of nearly half an electron. 

In addition we studied the effect of elongation of the CO distance on the top-position 
adsorption site. As a result we found only a very small elongation of 0.01—0.02 a.u. compared 
with the free molecule. 



VI. CONCLUSION 

Using a cluster approach within a relativistic density functional approximation as de- 
scribed above, we get results for the adsorption of the CO molecule on the Pt(lll) surface 
which are reasonably good. This statement in itself is very remarkable although we must 
admit that we are not able to show the convergence of these results with cluster size here. 
For one of the positions we have been able to include more Pt atoms in the cluster and 
in fact, the binding energy did not change any more. As expected, the binding energy is 
slightly higher than the experimental value of Yeo et al. j20| . 



Finally we note that the used method is able to describe the adsorption process. We are 
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able to calculate this important system with a quasi ab initio method (the only experimental 
result used in the calculation is the structure of the bulk) in a fully relativistic framework. 
In our calculations all relativistic d orbitals of Pt are treated dynamically. 
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Table Capture 

Table J). Results for the CO molecule: binding energy (BE), binding distance (dc-o), ion- 
isation potential (IP), frequency (ojc-o) of the molecular stretch vibration and dipole 
moment (/i). The dipole moment was calculated at the distance of the energy mini- 
mum. Assuming the experimental binding distance a value of /i = —0.052 a.u. was 
obtained. 

Table [□]. Results for the binding energy of CO on Pt(lll) in eV for the four adsorption 
sites. 

Table |T|. Results for the dist ance of the carbon atom to the surface (in a.u.) for the four 
bond positions. 

Table [IV]. Comparison of the results of the distance of the carbon atom to the next surface 
platinum atom (in a.u.). 

Table [V]. Comparison of theoretical results for the top-position. 

Table [VI| . Comparison of theoretical results for the hollow-position. 

Table |V1J| . Electronic charge of the adsorbed molecule and the consisting atoms as a result 
of a Mulliken population analysis for CO at Pt(lll). 
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Funktional 


RLDA 


B88/P86 


Experiment 


BE (eV) 


14.18 


13.25 


11.09 m 


dc-o (A) 


1.148 


1.153 


1.128 m 


IP (eV) 


15.01 


15.11 


14.01 ||] 


io C -o{cm' 1 ) 


2281.11 


2176.22 


2169.81 H 


\x (a.u.) 


-0.027 


-0.020 


-0.046 @ 



TABLE I: 



Figure Capture 

Figure |J. Potential energy curve of CO on Pt(lll) at the top position. The platinum 
cluster consists of 13 atoms: Pti3(7, 6). Solid curve is obtained with B88 exchange 
and BP86 correlation potential and dotted curve is obtained with PW91 exchange and 
correlation potential. 

Figure Potential energy curve of CO on Pt(lll) at the bridge position. The platinum 
cluster onsists of 7 atoms: Pty(4, 3). Solid curve is obtained with B88 exchange and 
BP86 correlation potential and dotted curve is obtained with PW91 exchange and 
correlation potential. 

Figure |3|. Potential energy curve of CO on Pt(lll) at the fee hollow position. The Platinum 
cluster consists of 12 atoms: Pti2(6, 6). Solid curve is obtained with B88 exchange 
and BP86 correlation potential and dotted curve is obtained with PW91 exchange and 
correlation potential. 

Figure |4j. Potential energy curve of CO on Pt(lll) at the hep hollow position. The plat- 
inum cluster consists of 13 atoms: Pti3(6, 7). Solid curve is obtained with B88 exchange 
and BP86 correlation potential and dotted curve is obtained with PW91 exchange and 
correlation potential. 
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Energy 
functional 


top- 
position 


bridge- 
position 


hep hollow- 
position 


fee hollow- 
position 


RLDA 
B88/P86 

Experiment [20] 


3.25 
2.30 

1.86 ±0.20 


3.15 
1.87 


3.35 
1.65 


3.25 
1.76 



TABLE II: 



Energy 


top- 


bridge- 


hep hollow- 


fee hollow- 


functional 


position 


position 


position 


position 


RLDA 


3.50 


2.60 


2.55 


2.55 


B88/P86 


3.57 


2.70 


2.65 


2.65 



TABLE III: 



Energy 


top- 


bridge- 


hep hollow- 


fee hollow- 


functional 


position 


position 


position 


position 


RLDA 


3.50 


3.69 


3.95 


3.95 


B88/P86 


3.57 


3.76 


4.03 


4.03 


Experiment 21 


, 22| 


3.50 ±0.19 


3.93 ±0.13 







TABLE IV: 



method 


v Pt -c[a.u.} 


E 6 (eV) 


Ref. 


Explanation 


RLDA 


3.50 


3.25 


this work 


13 Pt atoms, 2 layers (7,6) 


RGGA 


3.57 


2.30 


this work 


13 Pt atoms, 2 layers (7,6) 


CASSCF 


3.84 


1.44 


m 


Complete Active Space 


PP-Slab 


3.55 


1.45 




6 layers slab 


ZORA-DFT GGA 


3.59 


1.41 


m 


spin-orbit, 2 layers slab 


exp. 


3.49 ±0.19 




[21, 22] 





TABLE V: 
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method 


T Pt _c[a-u.} 


E b (eV) 


Ref. 


Explanation 


RLDA 
RGGA 
ORPP 

ZORA-DFT GGA 
exp. 


3.95 
4.02 
3.96 
4.12 
3.93 ±0.13 


3.25 
1.76 
2.39 
1.05 


tl 
tl 


lis 
lis 

[ 
[ 

21 


w 
26 
35 

5 | 


'ork 
r ork 

1 
1 

2j 


9 Pt atoms; 2 layers (6,3) 
9 Pt atoms; 2 layers (6,3) 
optimized relativistic pseudo potential 
pin-orbit, 2 layers slab 



TABLE VI: 





top- 


bridge- 


hep hollow- 


fee hollow- 


free 




position 


position 


position 


position 


molecule 


Oxygen 


8.375 


8.370 


8.359 


8.359 


8.348 


Carbon 


5.719 


6.072 


6.094 


6.075 


5.652 


Sum 


14.094 


14.442 


14.453 


14.434 


14.000 



TABLE VII: 




3.2 3.4 3.6 3.8 4 4.2 

distance C-surface in a.u. 

FIG. 1: 
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-1.5 




" 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 



distance C-surface in a.u. 
FIG. 2: 




" 2.2 2.4 2.6 2.8 3.0 3.2 3.4 



distance C-surface in a.u. 
FIG. 3: 



21 




' 2.2 2.4 2.6 2.8 3 3.2 3.4 



distance C-surface in a.u. 



FIG. 4: 
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